Light-front wavefunction dependence of the quark recombination 
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Abstract 

We present an extension of the recombination formalism to analyze the effects from the variation 
of the hadron wavefunctions. The hadron spectra are sensitive to the shape of the wavefunctions. 
However, when we fit the wavefunction parameters to the physical observables, such as the average 
charge radius, the final spectra are very similar each other. We discuss our numerical results in 

in 

' comparison with the published PHENIX and STAR data at RHIC. In the hadron spectra, the 

^ • recombination of thermal partons dominates at intermediate transverse momentum (Pt = 2 ~ 5 

Qh. GeV), and the fragmentation dominates at high P^ (> 5 GeV). The yield ratios and the nuclear 



modification factors for various hadron species are also estimated and compared to the experimental 
data. We present a new prediction on p/p and K~ /K^ ratios, including the jet quenching effects 
to the fragmentation mechanism. 
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I. INTRODUCTION 



Although no one doubt the existence of quarks and gluons, they have not yet been 
detected individually at the zero temperature. The current vigorous efforts at the Relativistic 
Heavy Ion Collider (RHIC) and the future plans of the Large Hadron Collider (LHC) may 
reveal the temperature dependence of the confinement mechanism and the chrial symmetry 
restoration Q]. The high-energy nuclear collisions compress and heat the heavy nuclei so 
much that their individual protons and neutrons may overlap and, in addition, a lot of 
pions may arise to ultimately create the quark-gluon plasma (QGP). The QGP may have 
existed ten millionths of second after the big bang and created the primordial matter of the 
universe. The RHIC and the future LHC may yield the QGP in the laboratory. It has been 
reported that the four experiments at RHIC already obtained distinguished results from the 
lower energy heavy ion collisions at CERN SPS Q, S, fl Q. The future LHC experiments 
(ALICE as well as CMS and ATLAS) would require theoretical predictions at the 30-fold 
energy increase from the RHIC. 

Among many others, the effects from the jet quenching and the bulk hadronization may 
be regarded as the important new results. Especially, the elliptic flow analysis revealed 
that the differential second-harmonic Fourier moment (^2) of the azimuthal distribution 
with respect to the reaction plane had a remarkable saturation property in the intermediate 
transverse momentum (Pt) range between 2 and 6 GeV for all hadrons including multi- 
strange baryons. This saturation effect and eventual decrease of V2 at high Pt have been 
qualitatively interpreted to be the results of partonic energy loss in an opaque parton system 
created by nuclear collisions P,!?!. Furthermore, the estimated V2 parameter as a function of 
Pt are scaled by the number of constituent quarks of particles. Together with an enhanced 
proton production in the intermediate Pt region, this agreement has been taken seriously 
as one piece of evidence for the quark recombination process and the presence of partonic 
collectivity at the early stage of a collision 0,0, [ifl, U, 12, 3]. 

In this work, we utilize the previous recombination formalism, and extend it to analyze the 
ight-front (LF) wavefunction dependence in the theoretical predictions from this formalism 
8|. Typical forms of the LF wavefunctions such as the Gaussian form and the power-law 
form are applied to this extended formulation. The numerical results are contrasted to each 
other, and compared with the available single invariant spectra by PHENIX and STAR 
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for various mesons and baryons. We also discuss the production ratios of various hadrons, 
including p/p and K~ jK^ ^ in the fragmentation region. While we include the jet quenching 
effects as others do, we get a rather distinguished results compared to the previous ones. For 
the high Pj^ regions, we get a dramatic suppression of the antiparticles p and K~ compare 
to the corresponding particles p and ^ respectively. 

The paper is organized as follows. In Sec. m we present the recombination formalism 
which is extended from the previous one to explicitly include the intrinsic transverse mo- 
menta of the constituents inside the hadron. Rather than an extensive review of the previous 
formalism, we focus on what has been extended from the previous model. In Sec. IIIH we 
present the numerical results of the Pt spectra for various mesons and baryons to contrast 
the results between the Gaussian form and the power-law form. The results are compared 
with the available experimental data from PHENIX and STAR collaborations. We also 
discuss the production ratios of various hadrons and the nuclear modification factor Rcp in 
this section. Conclusions and discussions follow in Sec. IIVI 

II. FORMULATION 

A. Recombination and Light-Front Wavefunction Dependence 

The current data from the RHIC experiments seem to indicate two distinguished mech- 
anisms of hadronization: (1) quark recombination for a rather low and intermediate Pt 
region and (2) quark fragmentation for a high P^ region. In this section, we present an 
extension of the recombination formalism to analyze the effects from the variation of the 
hadron wavefunctions. 

Introducing the density matrix p for the system of partons, the number of quark-antiquark 
states that one may interprete as mesons is given by 



where |M; P > is a meson state with the momentum P and the sum is over all combinations 
of quantum numbers such as flavor, helicity, and color of the valence quarks that contribute 
to the given meson M. Inserting complete sets of coordinates and using the Wigner function 
formalism, one can derive the formula for the invariant spectrum of the meson M as follows 




(1) 
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j^d^NM ^ fd^RP-u{R) f d^q 



dP^ 7s (27r)3 J (27r) 

d^RP-u{R) f dxP+d'^k 



C 



M 



(27r)3 J (27r)3 
xwaiR; xP^. k^)^Mix, k^)whiR; (1 - x)P+, -fc^), (2) 

where ^Yiil) — I d^r^Z{r,q) in the Wigner function formahsm and ^M{x,k±) = 
\iIjm{x, using the LF wavefunction of the meson ipAiix, k±). Here, x and k± are the mo- 
mentum fraction and the respective intrinsic transverse momentum of each quark. Similarly, 
the invariant spectrum of the baryon B can be obtained as follows ^: 

^d^Ne „ f d^RP ■ u{R) f dxiP+d%^ f dx2P^d%^ f dx^P+d%^ 



dP^ "7s (27r)3 7 (27r)3 J (27r)3 J (27r)3 

XWa{R] XiP+, ki^)Wb{R] X2P^, k2A_)Wc{R] X^P^ , k^^) 

X^B{xi,X2,X3,ki±,k2±,k3±), (3) 

where xi + X2 + X3 = 1 and fcij, + /i;2± + = 0. 

The previous work by Fries and collaborators used a factorized ansatz for the LF wave- 
function, for example [8|, 

^m{x, k±) = 0A/(x)fi(A;_L), (4) 

for mesons with a longitudinal distribution amplitude (f)M{x) and a transverse part Q{k±). 
However, we note that such factorization ansatz cannot be justified in free space since the 
LF wavefunction depends on the LF invariant mass of the particle, e.g., for the meson 
(m^ + kj_)/x + {ml + kjj/ll — x) (here the meson is composed of quark a and b), which 
cannot be factorized as Eq.Q. In general, the assumption of wavefunction factorization 
such as Eq.(jD) is not acceptable in free space because LF wavefunctions should be solutions 
of LF bound-state equations and the LF energy-momentum dispersion relation is rational, 
i.e. k~ = + m?)/k~^ for the particle with mass m. Both the LF kinetic energy (i.e. 
the LF invariant mass of the bound-state) and the LF kernel (or the inverse of the LF 
energy difference between the initial and intermediate states) involved in the LF bound-state 
equations are not factorizable due to the rational energy-momentum dispersion relations. 
Thus, the solutions of the LF bound-state equations cannot be factorizable and we do not 



integrate over k± in Eqs.© and Q but leave k± explicitly in the formulation. On the other 
hand, since it is not yet known if the LF bound-state solution in free space is also applicable 
to the dynamical recombination process in quark matter without any modification, we note 
that the factorization ansatz used in Ref. 10] may be equally valuable as one of the model 
calculations in this work. 

The usual parton spectrum at a given temperature is given by ^ 

where p and are the radial and the azimuthal angle coordinates, respectively. In addition, 
v{R) and f] represent the velocity four vector and the rapidity of the quark a, respectively. 
Here, 7a = exp(/ia/T) is the fugacity factor for each quark species a for w 
results from the statistical analysis for the hadron production at RHIC 
potential /i^'s are 9, 6.7, and —3.9 MeV for a = u(or d), s, and c, respectively. Note that 
Ref. [l^ gives resulting chemical potentials for isospin, strangeness and charmness as well 
as baryon chemical potential estimated by the statistical model at RHIC. Since statistical 
analysis of hadron production provides only chemical potentials of hadrons (not quarks), 
we obtained the quark fugacities by using the following formula for the fugacity of hadron i 



lich we adopt the 
141 ]: the chemical 



.1 



^^ = ^IsU^a' (6) 

a 

where is the number of quark species a in hadron i. In Eq.©, the fugacity is 
close to 1, as the isospin chemical potential and the assumed freeze-out temperature 
T are -0.96 MeV and 175 MeV, respectively |lj]. We assume that the temperature T for 
hadronization occurs at 175 MeV. The lattice QCD predicts that the phase transition occurs 
at T = 175 MeV at vanishing baryon chemical potential [16]. It should be reasonable that 
the temperature of the partonic phase is assumed to be the same as that of the phase 
transition. The space-time structure of the parton source in Eq.© is given by a transverse 
distribution /(p, (p) and a wide Gaussian rapidity distribution with a width A. Also, one 
may assume f{p,(f)) ~ ©(po ~ p) especially for the analysis of the central collisions. With 



^ The quark fugacity 7^ in this paper means Aq in Eq.(3) of Ref. [isf with the saturation factor 



one. 
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these assumptions, we find 



M 



Pt sinh rjT 
T 



X dx d'^k±\^ljMix,k±)\'^kMix,k±, Pt), (7) 
^0 Jo 

where V = tAt (t is the hadronization time and At is the transverse size) is the volume of 
the parton system and 

cosh r]T 



kM{x,k±,PT) = Ki 



T 



{^Jml + {xPt + k^y + J ml + {(1 - x)Pt - k^Y) 



(8) 



We note that the particular combination of Pt and fc_L for each constituent quark in Eq.(jHJ is 
consistent with the boost invariance of /cm in light-front dynamics. Extension to the baryon 
case is straightforward as 



dP^dy 



B 



V 

CbMt tt—ttt 27a7fe7c/o 



Pt sinh tjt 
T 



dxidx2 / d kij_d k2± 
Jo 



X \'iIjb{xi,X2, ki±, k2±)\'^kB{xi, X2, ki±, k2±, Pt), 



(9) 



and 



kB{Xl,X2,ki^,k2±,PT) 



+ A/m2 + {(1 - Xi - X2)Pt - (fclX + k2A_)Y}]. 



(10) 



In the following analysis, we take ^ as a free parameter to fit all invariant spectra of hadrons 
simultaneously for a given collision centrality. 

With this extension, we can now explicitly include the effect from the intrinsic transverse 
momentum of each quark, and var y th e form of the LF wavefunction such as the Gaus- 
sian form and the power- law form |l2|. In this analysis, we use the following typical LF 
wavefunctions for mesons and contrast the results between the two: 



'4'Gaussian^X, k±) — CXp 

and 

'4^Power—lawix, k±) 1/ 



X 



1 — X 



(11) 



— H 



X 



1 — X 



(12) 
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where /3, a, and n are the parameters that can be fixed from the physical observables such 
as the size and the mass spectrum of meson, etc.. The extension of Eqs.()lip and (|T^ to 
baryons is rather straightforward: 



1pGaussian{Xl,X2, ki±, k2±) = CXp 

and 

'ipPower-law{Xl,X2,kij_,k2±) = 1/ 



Xi 



X2 



+ 



Xi - X2 



(13) 



ml + fc^_L , ml + k1 



+ 



• (14) 



Xi X2 I - Xi- X2 

In this calculation, we used 260 MeV for the masses of u and d quarks and 460 MeV 
for the mass of s quark. In a relativized quark model with chromodyamics, spectra of 
both mesons and baryons have been well analyzed. As shown in typical references ([isl 
for mesons and 19] for baryons), the potentials among constituents such as confinement, 
hyperfine, spin-orbit, etc. work together to reproduce the hadron spectra comparable to the 
experimental values. For instance, proton mass was predicted as 960 MeV and the mass 
difference between nucleon and delta was obtained around 300 MeV while the light quark 
mass was taken as 220 MeV. The same light quark mass was used to predict the meson 
spectra which were overall in good agreement with data. These support our light-front 
quark model calculations (see e.g. Although we took the light quark mass 260 MeV 

as used in Ref. ^ for the present analysis, the essential predictions from a relativized quark 
model (or light-front quark model) remain intact. 

Just to illustrate how typical LF wavefunctions look like, we plot '>pGaussian{,x,kx.) for 
different (3"^ values in Fig. ^ As expected, the LF wavefunctions are symmetric around x = 
0.5, if the masses of constituent quark and antiquark are the same for pions. As increases, 
i'Gaussianix, k±) bccomcs broader in x as well as in kj_. When the mass of constituent quark 
and antiquark are not the same like K and D, ipcaussianix, fc_L)'s are clearly skewed in x. 

In order to constrain P, a, and n in Eqs.(fTT|) - (HH), the average values of k± are fixed 
by the measured average charge radius square < > of each hadron. If the experimental 
data for < > are not available, we adopt the calculated ones by a relativistic quark model 
2l[ [2^ . As an example, < > = 0.44 fm^ for pions, and the corresponding is 0.825 
GeV^ for ipGaussian- The average values of the charge radius square and the corresponding 
values of are summarized in Table HI In addition, the deduced and n of ippower-iaw are 
0.5 (1.53) GeV^ and 2 (6), respectively, for pions (protons). The left two panels of Fig. |21 



show the comparison of ipGaussian and ijjpower-iaw of pions by using the adjusted wavefunction 
parameters. They demonstrate that the LF wavefunctions are very similar in shape when 
the parameters are determined by some physical observables such as the charge radii of 
hadrons. However, if we use some arbitrary values for those parameters, the shape of LF 
wavefunctions can be quite different, which is demonstrated in the lower right panel of Fig.|21 
The importance of the proper choice of the LF wavefunction parameters in the hadron 
spectra are explained in Fig. El for vr's and protons. The invariant yields of the recombined 
hadrons are quite different for different sets of parameters. However, the hadron yields 
from the recombination process are quite similar for ipcaussian and i'power-iaw, once the 
wavefunction parameters are fixed by some physical observables (See solid vs. dashed lines 
in Fig. El). In the following analysis, we use only the Gaussian wavefunction with proper 



for each hadron. For the comparison, the wavefunctions used in Ref. 



are also included 



in Fig. El The wavefunctions used in the present analysis (solid lines) are lower than the 
factorized wavefunctions used in Ref. ^ shown by dotted lines especially for relatively low 
Pt region. Even in the logarithmic scale, the differences are as visible as the ones with 
arbitrary wavefunction parameters (dash-dotted lines). 

Since the essence of this work is to study the effect of the proper treatment of the LF 
wavefunctions to the recombination yields, we also compare the invariant spectra estimated 
by using full LF wavefunctions with those estimated by factorized wavefunctions. For the 
approximation that uses factorized wavefunctions for mesons, we tested ipGaussiani^, k±) of 
Eq.(fTT|) with < kl >= 0.088 GeV^, multiplied by fi(A?x) = exp (-A?!//?^) with p'^ = 0.176 
GeV^, which gave us the right < k'j_ > value. The resulting invariant spectrum for the 
recombined vr^'s is very similar to the PLl option in Fig. El Since it duplicates the result of 
PLl option, we do not display it explicitly in Fig. El However, this result indicates that the 
above factorized wavefunction may be as useful as the full LF wavefunctions developed in 
this paper for the phenomenological analysis of mesons, when the parameters are properly 
chosen. 

Similarly, for baryons, we tested i'Gaussian(yXi,X2,ki±,k2±) of Eq. (fT!?|l with < A;.fj_ >= 
0.0512 GeV2, muhiplied by n{k^) = exp [-{kf^ + kl^)/2(3^] with (3^ = 0.495 GeV^. How- 
ever, we could not find proper value for the factorized wavefunction to get the correct 
< k'j_ > = 0.0512 GeV^ since the results were not at all stable (too large standard devia- 
tions). Thus, we instead varied value in a wide range starting from 0.03 GeV^ all the way 
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even above 1 GeV^, and compared the resulting proton spectra. Results for P"^ below 0.03 
GeV^ couldn't be obtained due to a numerical instability. The results for above 1 GeV^ 
were about the same as the result of = 1 GeV^ and very stable as expected from the form 
of the above factorized wavefunction. Also, what we obtained for values from 0.03 GeV^ 
to 1 GeV^ was that the proton spectrum results were fairly insensitive to the (3"^ value. The 
resulting invariant spectrum of the recombined protons are shown by the coarsely-dotted 
line (typically for = 0.495 GeV^ but the result doesn't change much for other jS"^ values) 
in the bottom panel of Fig. El The recombined proton yield by the factorized wavefunction is 
significantly lower than the one by nonfactorizable wavefunction. This comparison indicates 
that it may be more significant in the baryon case than in the meson case what form of 
the wavefunction is taken for the prediction of the recombination process. Thus, the proper 
treatment of the wavefunction seems particularly important for the baryon production in 
the recombination process. 

We note that there are also other approaches for the recombination process: Refs. Q| 
and 23] considered soft-hard recombination with a covariant coalescence model and the 
fragmentation as a part of recombination, respectively. In addition, a recent work by A. 



Majumder, E. Wang, and X.-N. Wang 2J] discussed a derivation of the recombination model 
from field theory description of jet fragmentation. They noted that an ad hoc formulation 
of the recombination model is only valid under some strict conditions on the hadron wave- 
function. Therefore, we note that the present development is not the first one to consider 
the dependence of hadron wavefunctions in the recombination process. 

B. Fragmentation and Jet Quenching 

Inclusive hadron production by fragmentation at large momentum transfer can be de- 
scribed well by perturbative quantum chromodynamics (pQCD). In the framework of pQCD, 
the invariant yield of hadron h with momentum P is given by [8j 

where the sum runs over all constituent quark species a in h. For the spectrum of parton a 
with momentum pa = P/z at midrapidity, we use the parameterization by pQCD: 
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where the parameters C, B, and k, are taken from a leading order pQCD calculations j25|. 
K = 1.5 is taken in order to consider higher order corrections approximately ^. Note that 
the number of partons in different collision centralities are obtained by scaling the number of 
binary nucleon-nucleon collisions {N^u) or, equivalently, by the nuclear thickness function 
{Taa)- The probability that parton a fragments into hadron h is taken into account by the 
fragmentation function 

Da^h{z) = Nz\l - zY, (17) 

where the numerical values of A^, 7, and 5 are taken from the parameteriza- 
tion by Kniehl, Kramer, and Potter (Tabel 2 of Ref. and the website 
" |http : //www . desy . de/~poette r/kkp . f ' ) for fragmentation of pions, kaons, protons and 
antiprotons. We call this parametrization as the KKP parametrization. A fragmentation 
'unction is taken from the work by de Florian, Stratmann and Vogelsang (Tabel 1 of Ref. 



Finally, the energy loss of energetic partons (so called iet quenching), especially, in central 
collisions is considered by the following parameterization 0, 

ApT(&,PT) = e(6)Vp?^^, (18) 

where Ra is the radius of nucleus A, {L){b) is the geometrical factor of the overlap zone of 
two nuclei, and e{b) is the energy loss parameter of the hot medium with impact parameter 
b. The detailed functional forms of {L){b) and e(6) are the same as Ref. j^: 

and 



but, practically speaking, it is reasonable to assume that {L){b) ^ Ra and e{b) ~ eo = 0.82 
GeV^/^ for the most central collisions as 6 — > 0. 



III. RESULTS 



A. Invariant Spectra 



Fig.|3]shows the numerical results on the invariant spectra of various mesons at midrapid- 
ity for central Au + Au collisions at ^/JNN = 200 GeV. In Fig.lU we compare our calculations 
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for the meson spectra with the pubhshed PHENIX and STAR data [30|, ISjJ, |32, |33|, |3J| up to 

Pt = 10 GeV in order to show the overall shapes, especially, the transition regions near 5 
GeV. The neutral pion spectrum was measured by PHENIX up to 10 GeV in Pt, but lacks 
data in a low Pt region. In contrary, the charged pions were measured only at low Pt up 
to 3 GeV with high precision. However, the high Pt spectra of charged pions are expected 
to be very similar to those of neutral pions. For charged kaons, PHENIX measured up to 
2 GeV, and STAR measured up to about 0.7 GeV in Pt- In general, the data by PHENIX 
and STAR agree quite well in the overlapped phase space except 0: the STAR data is about 
a factor of three larger than the PHENIX data. 

In Fig. m the dashed and dotted lines represent the model calculations from the re- 
combination and the fragmentation, respectively, and the solid lines are the sum of two 
contributions. In the 7r° spectra, the two distinguished Pt regions of hadron production 
are manifest. Although the transient Pt depends on the particle species, the recombination 
process is dominant between ~ 2 and 5 GeV, and the fragmentation is dominant above 5 
GeV. Our calculation is not expected to reproduce the hadron spectra below about 2 GeV 
in Pt- In such a very low Pt region, the calculation underestimates the experimental data 
significantly, implying that other processes like the transverse flow, the secondary decay of 
hadron resonances, and the binding energy effect become important. Our calculation re- 
produces the measured meson Pt spectra larger than 2 GeV reasonable well, including the 
strange mesons. Note that we do not plot the fragmentation contribution for due to the 
lack of the fragmentation function. 

Fig. Elshows the numerical results on the invariant spectra of baryons up to Pt = 10 GeV 
at midrapidity for central Au + Au collisions at ^/snn = 200 GeV. Fig. El also compares 
;he calculations with the published experimental data with the open circles by PHENIX 



3l[ Issj l and the solid triangles by STAR js^, 36, Is^. The left most column of Fig. |S1 is 



for protons and antiprotons; PHENIX and STAR measured up to about 5 and 1.2 GeV, 
respectively. Note that the published Pt spectra for protons and antiprotons by STAR 
j^d] are about 40 % higher than those by PHENIX Slj- This difference comes from the 
fact that the contributions from the A and S'' decays are removed only for PHENIX. For a 
fair comparison between the two data sets, the p and p spectra by the STAR collaboration 
are scaled by 0.6 in Fig. After scaling down, the STAR spectra agree quite well with 
the PHENIX spectra in the overlapped phase space. The present model reproduces the 
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measured proton and antiproton spectra reasonably well. The model also predicts that the 
transient Pt for baryons from recombination to fragmentation is somewhat higher than that 
for mesons. 

In Fig. 13 the experimental invariant spectra of A + 2", H^, Vt^ , and their antiparticles are 
for y^s^N = 130 GeV. But all model calculations are for ^/SNN = 200 GeV because all input 
parameters of the model calculations are available only for ^/snn = 200 GeV. Furthermore, 
due to the lack of the fragmentation functions, we do not plot the fragmentation contribution 
for S and Q. Because of the difference in beam energy, the model overestimates the yields 
of A + S", , and their antiparticles, and the discrepancy is larger for a larger number 
of strange quarks in a given baryon. 



B. Yield Ratios 



One of the most interesting data from RHIC is the yield ratio of protons(or antiprotons) 
to pions at the intermediate transverse momentum region (2 < < 5 GeV) in central 
heavy-ion collisions. The p/ir and p/ir ratios rise steeply with Px up to about 2.5 GeV, but 
levels off at about 1 and 0.7, respectively, in 2.5 < Pr < 5 GeV for the most central 10 % Au 



+ Au collisions 



31 



|. At Pt > 2 GeV, p/ir and p/ir for peripheral collisions are similar 
to those for elementary p + p and e^e~ collisions, and the ratios increase from peripheral 
to central collisions. Since the hydrodynamic model, which had been rather successful in 
describing the low Pt hadron spectra, could not explain the centrality dependence of limiting 
values, a recombination mechanism of hadronization at intermediate Pt was proposed as a 
possible resolution 10, 23, The recombination process naturally explains that 



the p/tt and p/n ratios level off in 2 < Pp < 5 GeV, and fall sharply near P^ ~ 5 GeV 
where the fragmentation takes over the recombination. Similar trends can also be found in 
the present calculation. The top row of Fig. IHl shows the results from our calculations for 
the p/vr", p/tc^, and p/p ratios in comparison with the published PHENIX data [31]. As Pt 
increases, the p/ir^ and p/ir^ ratios rises, reach the maximum values around 3 GeV, decrease 
sharply, and, finally, become constant at about 0.1 for Py > 6 GeV. In addition, the p/p 
ratio is almost constant at about 0.9 for P^ < 5 GeV. However, it also decreases with Pt, 
and become almost constant at about 0.1 for Pt > 7 GeV, which is very different from the 
previous calculation by Fries et al. (the dashed line in Fig. ^ P] . Although there is no 



12 



dispute on KKP parametrizations [2y| for the gluon fragmentation, the quark fragmentation 
is a problem because the KKP fragmentation functions are not fully flavor seperated and 
one has to make additional assumptions to seperate contributions from different flavors. 
Incidentally, the recent STAR data on identifled hadrons 41] reveal a poor job of the KKP 
fragmentation functions for p and p yields, possibly due to the lack of flavor seperation. 
It seems that the sea quark contributions in and p/p are particularly problematic. 
For instance, as stated in Ref. 2^, the d quark in does not behave sea-like contrary 
to expectations. Thus, even a slight difference in handling sea quark contributions could 
make a large difference in the predictions of ratios for K^/K^ and p/p. Our results are 
based on maintaining expected smallness of sea quark contributions consistently not only 
in the pion case but also in other hadron cases. Our low ratio for Pt > 7 GeV in p/p 
of Fig. El is consequently due to the dominance of valence contributions. It is a fact that 
incident heavy ions possess valence quarks, but not antiquarks. In other words, the charge 
conjugation symmetry is already broken in RHIC environment due to the initial nuclei carry 
only nucleons (not antinucleons). Since the baryon number (or, equivalently, the quark 
number) must be conserved throughout the reactions, more protons than antiprotons are 
expected in the fragmentation region. 

For more comparisons on the hadron yield ratios, the bottom row of Fig. [Blshows /tt^, 
K~/tt~, and K~/K^. Although the measured Pt range of is limited, the present 
estimates are in reasonable agreements with the data. The K^/tt^ and K^/ir^ ratios 
increase with Pt, and reach maximum around Pt = 3 GeV. If Pt further increases, /-k^ 
and /t\^ decrease, and level off at some constants. The Pt dependence of the / 
ratio is very similar to that of p/p. Especially, we note that the present results on /-r' 
and K~ / at high Pt region, where the fragmentation is dominant, are quite different 
from the previous model calculations by Fries et al. [8]. As discussed above, even a slight 
difference in handling sea quark contributions could make a large difference in the predictions 
of ratios for K~ / as well as /-k^ . The forthcoming RHIC data at high Pt, e.g., the 
PHENIX data with newly installed aerogel detector, ma y help to further clarify the flavor 
seperation issue in the KKP fragmentation functions 
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C. Nuclear Modification Factors 



Another important feature of the RHIC data can be identified in the nuclear modification 
factor Rcp, which is defined by the Ncou scaled central to peripheral yield ratios: 

'W' ■) Jcentral I ^ t\j central ^ 

The RHIC experiments observed that the Rcp parameters of various mesons in > 2 GeV 
in central collisions were suppressed with respect to the Ncou scaled p + p and peripheral 
collision data. Moreover, the suppression in the intermediate transverse momentum region 
(2 < < 4 GeV) was only for mesons, but not for baryons. The experimental Rcp param- 
eter of protons in intermediate transverse momentum region is unity, which is completely 
consistent with N^oii scaling. The Rcp of A and A are also close to unity in an intermediate 
Pt region, but somewhat smaller than protons. 

One of possible explanations for the suppression of hadron yields at high Pt and a dis- 
tinguished behavior of mesons and baryons at the intermediate Pt region is the combined 
effect of recombination and fragmentation. Fig.[7|shows the estimated Rcp parameters of vr, 
p, and A + A as a function of Pt- For pions, we plot charged and neutral pions together, as 
almost no difference is expected from the present model. In the present model calculation, 
we scaled the hadron yields due to the fragmentation by the number of binary collisions. 
For the recombination part, the ratio of ^]^7a in Eq.((Zj) for peripheral collisions to that in 

a 

central collisions was assumed as 40 % of the number of participant {Npart) ratio: 

/^TT \ peripheral 

\ J_J_ ji^peripheral 

^1 ATcentra.l ' (22) 



central TUcentral ' 

^ part 



where Ci = 0.4 fitted by the Rcp parameters of vr, p, and A + A, simultaneously, by fixing 
the temperature at 175 MeV as it is almost independent of the collision centrality 
Since the quark fugacities are also almost constant except for the very peripheral collisions 



42j |. the factor ci mostly reflects the effect of volume. As a result, the fact that c < 1 



is understandable as the flow velocity is larger for more central collisions. The agreement 
between the present calculations and the experimental data are reasonable for all considered 
hadron species. 
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IV. CONCLUSIONS 



We have presented an extended formalism of the recombination model to analyze the 
effects from the variation of the hadron's hght-front wavefunctions. Two different functional 
forms of the light-front wavef unction, which are the Gaussinan form and the power-law 
form, are tested in detail. The hadron spectra are indeed sensitive to the shape of the wave- 
functions. However, when we fit the wavcfunction parameters to the physical observables, 
such as the average charge radius, the final spectra are very similar each other. We discuss 
our numerical results in comparison with the published RHIC data, especially, from the 
PHENIX and STAR collaborations. In the hadron spectra, the recombination of thermal 
partons dominates at the intermediate transverse momentum region between 2 and 5 GeV, 
and the fragmentation dominates at high Pt larger than 5 GeV. The yield ratios and the 
nuclear modification factors for various hadron species are also estimated. In general, the 
present model, which combines the recombination and fragmentation processes, are quite 
consistent with the experimental data. We have also discussed new predictions an. p/p and 
K~ / K'^ ratios, including the jet quenching effects to the fragmentation mechanism. 
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FIG. 1: Shapes of the Gaussian wavefunctions as functions of x and k_\_ for different (normaUza- 
tion is not performed in this figure). The top, middle, and bottom rows represent the wavefunctions 
for vr, K, and D, respectively. 
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FIG. 2: Comparison of the pion wavefunctions between the Gaussian form and the power-law 
form. The Gaussian wavefunction with 0^ = 0.825 GeV^ (top) and the second order power-law 
wavefunction with = 0.5 GeV^ (bottom left) are adjusted to the average charge radius of pions. 
But the first order power-law wavefunction with = 0.825 GeV^ (bottom right) is not adjusted. 
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FIG. 3: Comparison of the invariant spectra of tt~^ (top) and p (bottom) by the recombination 
process for various assumptions on the wavefunction parameters. Left panels show the spectra 
in hnear scale in order to emphasize the difference in relatively low Pt region, while right panels 
show them in log scale for the comparison of the overall shapes. G represents the Gaussian 
wavefunction used in this paper. PLl represents the power-law wavefunction whose parameters are 
adjusted by the known average charge radius, whereas PL2 represents the power-law wavefunction 
with arbitrary values for the parameters. The coarsely-dotted lines are calculated by a factorized 
Gaussian form of the LF wavefunctions with = 0.495 GeV^ (see text for details). For the 



comparison, the wavefunctions used in Ref. 



are also plotted by the dotted lines. 
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FIG. 4: Invariant spectra of mesons at midrapidity for central Au + Au collisions at y^F/vTV = 200 
GeV. The dashed and dotted lines represent the model calculations from the recombination and 
the fragmentation, respectively. The solid lines are the sum of two contributions. The open circles 



are the published data by the PHENIX collaboration 



by the STAR collaboration 



2L, 



and the solid triangles are those 
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FIG. 5: Invariant spectra of baryons at midrapidity for central Au + Au collisions at ^/sJfN = 
200 GeV. The dashed and dotted lines represent the model calculations from the recombination 
and the fragmentation, respectively. The solid lines are the sum of two contributions. The open 



circles are the published data by the PHENIX collaboration 

"3,137 



those by the STAR collaboration 



33 



31 



35l | whereas the solid triangles are 



. For a fair comparison between two sets of the data, 



the published p and p spectra by the STAR collaboration [33| are scaled by 0.6, which removes 
the contribution by the weak decays of A, T,^, and their antiparticles. Note that the experimental 
invariant spectra of hyperons are for ^/sNN = 130 GeV, whereas all model calculations are for 
y^i]^ = 200 GeV. 
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FIG. 6: Calculated hadron yield ratios (solid lines) as a function of Pt in comparison with the 
PHENIXdata 3- 

For the comparison, we also show the model calculations by Fries et al., in 
p/p, K~ /tt~, and /K^ by dashed lines 
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FIG. 7: Nuclear modification factors Rcp as a function of Pt for tt, p, and A + A. The solid lines 
are from the present calculation. 
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